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ABSTRACT 

bJQ' 

3 : 

' Two possible explanations for the low luminosity of the supermassive black hole at 

CO ' the center of our galaxy are (1) an accretion rate of order the canonical Bondi value 

(~ 10~ 5 Mq yr -1 ), but a very low radiative efficiency for the accreting gas or (2) an 
fNj ■ accretion rate much less than the Bondi rate. Both models can explain the broad-band 

spectrum of the Galactic Center. We show that they can be distinguished using the 
linear polarization of synchrotron radiation. Accretion at the Bondi rate predicts no 
\ linear polarization at any frequency due to Faraday depolarization. Low accretion rate 

models, on the other hand, have much lower gas densities and magnetic field strengths 
close to the black hole; polarization may therefore be observable at high frequencies. If 
confirmed, a recent detection of linear polarization from Sgr A* at ^ 150 GHz argues 
for an accretion rate ~ 10~ 8 Mq yr , much less than the Bondi rate. This test can be 
q . applied to other low-luminosity galactic nuclei. 

a: 

Subject Headings: accretion, accretion disks — Galaxy: center — polarization 

1. Introduction 



In this paper we discuss the effect of Faraday depolarization on synchrotron radiation in spher- 
ical accretion flow models of low-luminosity galactic nuclei (see also Bower et al 1999ab; Agol 2000). 
We focus on the radio source Sgr A* at the Galactic Center, but our results can also be applied 
to other systems (see §4). This paper was motivated by a possible detection of linear polarization 
from Sgr A* (Aitken et al. 2000). 

The Bondi accretion rate onto the supermassive black hole at the center of our galaxy is 
estimated to be ~ 10~ 4 — 10~ 5 Mq yr -1 (e.g., Coker & Melia 1997; Quataert, Narayan, & Reid 
1999), implying a luminosity of ~ 10 41 ergs s _1 if the radiative efficiency is ~ 10%. This is roughly 
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5 orders of magnitudes larger than the observed luminosity (see Narayan et al. 1998 for a recent 
compilation). Comparable "discrepancies" are obtained for massive elliptical galaxies in nearby 
X-ray clusters (e.g., Fabian & Rees 1995; Di Matteo et al. 1999). 

One explanation for the low luminosity of nearby supermassive black holes is that they accrete 
via an advection-dominated accretion flow (ADAF), in which most of the dissipated turbulent 
energy is stored as thermal energy rather than being radiated (e.g., Rees et al. 1982; Narayan & 
Yi 1994, 1995; Abramowicz et al. 1995). In such models the accretion rate is of order the Bondi 
rate while the radiative efficiency is extremely small (~ 10~ 6 for Sgr A*). 

Another explanation for very low luminosity accreting systems is that the Bondi accretion 
rate estimate is inapplicable (e.g., Blandford & Begelman 1999; Gruzinov 1999). In particular, 
numerical simulations of non-radiating accretion flows with small values of the dimensionless vis- 
cosity parameter a find that the gas density scales with radius as p cx r -1 / 2 rather than the 
canonical Bondi/ ADAF scaling of p cx r -3 / 2 (Stone, Pringle, & Begelman 1999; Igumenshchev & 
Abramowicz 1999, 2000; Igumenshchev, Abramowicz, & Narayan 2000). Narayan, Igumenshchev, 
Sz Abramowicz (2000) and Quataert & Gruzinov (2000) explained these simulations in terms of 
a "convection-dominated accretion flow" (CDAF). In such a flow angular momentum is efficiently 
transported inwards by radial convection, nearly canceling the outward transport by magnetic 
fields. This strongly suppresses the accretion of matter onto the black hole. 

Broad band spectra have thus far had difficulty distinguishing between these explanations for 
the low luminosity of nearby supermassive black holes. For example, Quataert & Narayan (1999; 
hereafter QN) showed that accretion at much less than the Bondi rate could produce spectra quite 
similar to ADAF models. In this paper we show that linear polarization observations in the radio to 
sub-mm can provide a sensitive probe of the accretion rate onto the black hole, and help distinguish 
between degenerate spectral models. 

In the next section (§2) we present simple estimates of the physical parameters of the accretion 
flow relevant for our analysis. We then discuss Faraday depolarization in spherical accretion flow 
models of Sgr A* (§3). In §4 we compare these predictions with observational constraints on the 
linear polarization of Sgr A* and summarize our results. We also generalize our analysis to other 
low-luminosity galactic nuclei. 

Throughout this paper, we focus on accretion models of Sgr A*. An unresolved jet or outflow 
may, however, dominate the observed emission (e.g., Falcke, Mannheim, & Biermann 1993; Lo et 
al. 1998; Falcke 1999); this is briefly discussed in §4. 



2. Plasma Parameters for Sgr A* 
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Stellar kinematics show that there are pa 2.6 x 1O 6 M within pa 0.015 pc of the Galactic 
Center (Eckart & Genzel 1997, Ghez et al. 1998), centered on the radio source Sgr A* (Menten 
et al. 1997). The most plausible explanation is that Sgr A* is a ~ 2.6 x 1O 6 M0 accreting black 
hole. Sgr A* is believed to accrete the winds from nearby (~ 0.1 pc) massive stars (Krabbe et al. 
1991). The Bondi accretion rate of these winds onto the supermassive black hole is estimated to 
be pa 10" 4 - 10" 5 M yr" 1 (e.g., Coker & Melia 1997; Quataert, Narayan, & Reid 1999). 

If the accretion rate close to the black hole is of order the Bondi value the gas density near 
r ~ 1 is n pa 10 9 — 10 10 cm" 3 (since v r pa c near the horizon). The corresponding magnetic field 
strength, assuming rough equipartition with the nearly relativistic protons, is B pa 2 x 10 3 G. At 
such magnetic field strengths, relativistic electrons cooling by synchrotron radiation would have a 
cooling time much less than the inflow time of the gas. In order to not overproduce the observed 
radio to sub- mm luminosity of Sgr A*, the bulk of the electrons must therefore be marginally 
relativistic, with T e ~ 10 9 — 10 10 K. These plasma parameters (n, B, T e ) describe Bondi and ADAF 
models of Sgr A* (e.g., Melia 1992, 1994; Narayan, Yi, & Mahadevan 1995; Narayan et al. 1998). 
In such models the electrons are assumed to be adiabatically compressed from large radii in the 
accretion flow, with virtually no additional turbulent heating. 

QN showed that accretion at much less than the Bondi rate could also produce the observed 
high frequency emission from Sgr A*, provided the electrons were much hotter than in standard 
ADAF models (see their Table 2 and Fig. 8b). A simple explanation for this result can be obtained 
by applying the Burbidge (1958) estimate to Sgr A*. We consider synchrotron emission from a 
sphere of radius R containing relativistic electrons with a temperature kT e = ^m e c 2 . We take the 
electron heating rate to be comparable to the net turbulent (magnetic) heating rate. As can be 
confirmed a posteriori, the synchrotron cooling time is S> the inflow time of the gas. The electron 
energy density is then similar to the magnetic energy density 

n-fm e c 2 pa — . (1) 

The frequency of peak synchrotron emission and the synchrotron luminosity are given by 

p R 

J.-0.17 2 — (2) 

m e c 

and 

L pa a T cB 2 j 2 R 3 n, (3) 

where o~t is the Thomson cross section. 

We express n, 7, and B in terms of R, u, and L (see also Falcke 1996; Beckert & Duschl 1997) 

7 »3.2(^^y /7 » 100, (4) 

n pa _f Pa 10 6 cm" 3 , (5) 
7 5 A 2 r e v ' 
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and 

B w y / 8^m^yn pa 45 G, (6) 

where A = c/V is the wavelength and r e = e 2 /(m e c 2 ) is the classical electron radius. For the 
numerical estimates in equations we have used the observed values for Sgr A*. The peak 

synchrotron frequency is at v sa 10 3 GHz with a luminosity of L ~ 10 36 ergs s _1 (e.g., Serabyn et 
al. 1997). In spherical accretion models, this high frequency emission arises from very close to the 
black hole, so we have taken R ~ R g ps 10 12 cm. 

Equation (|l~5|) gives a density close to the black hole of n ~ 10 6 cm -3 ; the implied accretion 
rate is then ~ 1O _8 M0 yr _1 , three to four orders of magnitude smaller than the Bondi value. 

The thermal blackbody emission at frequency v from a sphere of radius R is 

L t = 2vri/ 3 7m e 47r J R 2 sa 10 37 ergs s _1 , (7) 

where the numerical estimate is for our fiducial parameters. This comparison shows that the 
synchrotron emission becomes optically thin below the peak frequency, near v ~ 300 GHz. At 
lower frequencies the emission is self-absorbed. 

The above considerations show that both low (~ 10~ 8 Mq yr _1 ) and high (~ 10 -5 — 10 -4 Mq yr -1 ) 
accretion rate models can explain the observed sub- mm "bump" in Sgr A*. Such models can be 
distinguished by comparing the observed brightness temperature and/or radio image as a function 
of frequency with the theoretical predictions (see, e.g., Ozel, Psaltis, & Narayan 2000). This test 
has been difficult to implement because interstellar scattering significantly broadens the image of 
Sgr A*.0 In the next section we show that the linear polarization of Sgr A* at high frequencies 
provides an additional discriminant. 



3. Faraday Depolarization 

The anisotropic index of refraction of a magnetized plasma leads to a frequency-dependent 
rotation in the position angle, 8, of linearly polarized electromagnetic waves, 

9 = RMX 2 , (8) 

where RM is the rotation measure. This can lead to significant depolarization of intrinsically 
linearly polarized synchrotron emission. 



2 Recent detections of Sgr A*'s intrinsic size (Lo et al. 1998; Krichbaum et al. 1998) still have sufficient uncertainties 
that a range of theoretical models are allowed. 
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For a "cold" non-relativistic plasma, RM is given by (e.g., Rybicki & Lightman 1979) 

RM = — / (fl-Bn 

2vrm 2 c 4 7 

/rcid 
cfl-Bn— (9) 
m 2 

where d\ is the differential path length from the observer to the source. In the Appendix we show 
that the rotation measure for an ultrarelativistic thermal plasma is given by 

RM^ = e -—, I dl • Bn b " 

2.63xl0- 13 x /fl-Bn^^ (10) 



7 27rm 2 c 4 J 2 7 2 

log 7 rad 
2 7 2 m 

where 7 = kT e jra e (? . A comparable expression is obtained for a power law distribution of rela- 
tivistic electrons, with 7 replaced by "fmin, the minimum Lorentz factor of the electrons (Jones & 
O'Dell 1977). In what follows, we define RM(r) to be the contribution to the net rotation measure 
from radii within dr ~ r of radius r in the accretion flow. 



3.1. ADAF/Bondi Models 

In spherical accretion flow models, higher frequency radio emission arises from closer to the 
black hole, where the electron temperature and magnetic field strengths are the largest; this is also 
true for jet models (e.g., Falcke 1999). Ozel et al. (2000) show that in ADAF models of Sgr A* 
the synchrotron emission at frequency v = 100 v\qq GHz arises from a radius r„ ~ 20 iq~oo 9 ( see 
their Fig. 5).Q This radius defines the r = 1 surface of the synchrotron emission. For smaller radii 
the emission is self-absorbed while for larger radii it is optically thin. Faraday rotation is only 
important for r ^ r u , where the photons "free stream" out of the accretion flow. 

In Bondi/ADAF models Faraday rotation is so strong that the synchrotron emission is com- 
pletely depolarized by the plasma within dr ~ r of r„, i.e., in the vicinity of the r = 1 surface 
where it is emitted (depolarization and emission are thus virtually co-spatial). Taking n oc r~ 3/2 
and B oc r -5 / 4 ; the rotation measure scales roughly as RM oc r -7 / 4 7 -2 . The relativistic suppres- 
sion of the rotation measure is small in all models which have an accretion rate comparable to the 
Bondi rate, because the electrons must then be at most marginally relativistic (§2). The rotation 
measure as a function of radius is thus given by (see also Bower et al. 1999ab) 

RM w 10 13 r~ 7/4 rad uT 2 . (11) 



3 Falcke (1999) finds a similar expression for r„ in the jet model. 
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In equation (11) we have assumed that the magnetic field is in rough equipartition with the gas 
pressure, has a significant component along the line of sight, and has a coherence length I ~ r; 
for I <C r, RM is reduced by ~ (^/r) 1 / 2 . The normalization in equation (11) is set by the Bondi 
accretion rate. 

The large RM in ADAF/Bondi models leads to a significant rotation in the position angle of 
linearly polarized waves. Photons of frequency v emitted at radius r v undergo Faraday rotation 
through an angle 

6 V « \ 2 RM(r v ) « lOVoK 7/4 ~ lOVoo 43 rad ' ( 12 ) 
where the last estimate uses Ozel et al.'s (2000) fit to r v {y). 

These rotation angles are so large that the synchrotron emission in ADAF/Bondi models of 
Sgr A* is completely depolarized by Faraday rotation. For example, in a simple uniform source 
model, the observed polarization is oc 0" 1 (Pacholczyk 1970). In general, the observed polarization 
depends on the rotation measure power spectrum, but is <C 1 for 6 V S> 1 (e.g., Tribble 1991). 



3.2. M < M B ondi 



If the accretion rate onto Sgr A* is much less than the Bondi rate, significant polarization may 
be observable at high frequencies; we show this using an order of magnitude estimate. 

For the M ~ 10~ 8 Mq yr -1 model of §2 the rotation measure calculated using equation ( |j~0|) 
is RM w 10 3 rad m 2 near r ~ 1. Moreover, if n oc r x / 2 , as in CDAF models, the magnetic field 
scales as B oc r~ 3 / 4 and 

RM w 10 3 r~ 1/4 (y^j rad m~ 2 . (13) 

The variation of the electron Lorentz factor with radius is somewhat uncertain, but we expect 
roughly 7 oc r _1 , so that the electrons become non-relativistic by r ~ 10 2 . Equation (|l~3| ) then 
shows that RM has its maximal value at large radii, r ~ 10 2 , where RM ~ 3 x 10 6 rad m~ 2 . 

Equation (|l~3|) demonstrates that, in contrast to ADAF models, there is no depolarization of 
synchrotron emission at small radii in models with accretion rates much less than the Bondi rate; 
RM is negligible in the region where the synchrotron emission is produced. Depolarization can still 
be important, however, because observed photons experience different Faraday rotation at large 
radii, r ^ 10 2 , on their way out of the accretion flow (e.g., Bower et al. 1999ab). 



4 One way of evading this conclusion is to posit that the magnetic field is sufficiently tangled (£ <§C r) to decrease 
RM to 10 6 rad m -2 . This tangling would, however, also eliminate linear polarization. 
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Spatial variation in the rotation measure will depolarize Sgr A* at frequencies for which 59 = 
X 2 5RM £ 7T, i.e., for 

^100 Ln. , -J GHz, (14) 



,10 6 rad m 2 , 

where 5RM is the difference in the rotation measure for photons of a given frequency which travel 
through different parts of the accretion flow.[] Quantitative calculations of depolarization by dif- 
ferential Faraday rotation are uncertain; two points are, however, clear: (1) At low frequencies, 
<C 100 GHz, Sgr A* is easily depolarized at r P> 10 2 . The required 5RM is <C 10 6 rad m -2 , orders 
of magnitudes smaller than the values of RM obtained at r ~ 10 2 — 10 4 . (2) Emission above ~ 100 
GHz can plausibly be linearly polarized if the accretion rate onto Sgr A* is much less than the 
Bondi rate. In particular, equations (13) and flH| ) show that for M <C MBondb emission above 
~ 100 GHz is not depolarized propagating out of the accretion flow. 



4. Discussion 



ADAF/Bondi models assume that the accretion rate onto Sgr A* is of order the Bondi rate 
(~ 10~ 4 — 1O~ 5 M0 yr _1 ) and that the radio to infrared emission is produced by synchrotron emission 
from marginally relativistic electrons (T e ~ 10 9 — 10 10 K). In such models the rotation measure 
is <; 10 10 rad m -2 inside ~ 100 Schwarzschild radii where the synchrotron emission is produced. 
ADAF/Bondi models thus predict that Sgr A* should be depolarized by Faraday rotation over the 
entire radio to infrared spectrum, and should have nearly zero linear polarization. 

The theoretical arguments summarized in §1 propose that the accretion rate onto Sgr A* is 
much less than the Bondi rate. We have described one such model, in which the electron heating 
rate is of order the rate of change of the magnetic energy density. For an accretion rate ~ 10 3 times 
smaller than the Bondi rate, i.e., ~ W~ 8 Mq yr _1 , and with relativistic electrons with 7 ~ 100, this 
model can explain the observed high frequency emission from Sgr A*. Moreover, it predicts that 
the rotation measure in the accretion flow is much smaller than in ADAF/Bondi models. This is 
because the gas density and magnetic field strength close to the black hole are much smaller, and 
because the electrons are relativistic (RM oc 7 _2 log7 for 7 ^> 1; see §3 and the Appendix). The 
maximal contribution to the rotation measure comes from ~ 10 2 — 10 3 Schwarzschild radii, where 
RM ~ 10 6 rad m~ 2 . 

Rotation measures of ~ 10 6 rad m -2 can depolarize Sgr A* at v <C 100 GHz by differential 
Faraday rotation; photons of a given frequency travel through different rotation measures on their 



5 To be precise, 8RM is the difference in the rotation measure at r ~ 100 — 10 4 on scales of r„, the source size 
at frequency v. This is difficult to calculate analytically, but could be determined from future MHD simulations of 
non-radiating accretion flows. 
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way out of the accretion flow. Following Bower et al. (1999ab), we believe that this accounts for 
the ^ 0.2% linear polarization of Sgr A* at low frequencies (from ~ 4 to ~ 23 GHz; see Bower et al. 
1999ab)|] it is less clear, however, that it can account for Bower et al.'s (1999b) limit of ^ 1% linear 
polarization at 86 GHz (see below). In fact, rotation measures of ss 10 6 rad m 2 are insufficient 
to depolarize emission above ~ 100 GHz. As a result, in models with accretion rates much less 
than the Bondi rate, <; 100 GHz emission is not depolarized propagating out of the accretion flow; 
intrinsically polarized synchrotron emission may therefore be observable at high frequencies. 

The above considerations show that the linear polarization of Sgr A* at high frequencies 
provides a means of distinguishing between accretion at the Bondi rate, and accretion at a much 
smaller rate. In fact, Aitken et al. (2000) report a possible detection of ~ 10% linear polarization 
from Sgr A* between 150 and 400 GHz. If confirmed, these observations require an accretion rate 
onto Sgr A* much less than the Bondi rate, roughly M ~ 10~ 8 AfQ yr _1 . 

One difficulty in interpreting Aitken et al's results is the large beam (~ 20") of the SCUBA 
camera on the JCMT. This large beam forced Aitken et al. to subtract out free-free and (polarized!) 
dust emission in order to isolate the flux and polarization of Sgr A*. Future high resolution 
polarimetry at mm wavelengths is clearly necessary to further address this important issue. 

Aitken et al. find that the position angle of Sgr A* changes by ^ 10° between A = 0.135 cm 
and A = 0.2 cm; at face value this implies RM ^ 10 5 rad m~ 2 , somewhat smaller than the values 
of ~ 10 6 rad m~ 2 in our model. This assumes, however, that the intrinsic position angle of Sgr A* 
is the same at A = 0.135 cm and A = 0.2 cm, which need not be the case. Moreover, our estimates 
of RM are actually upper limits, since they assume (1) equipartition magnetic fields aligned along 
the line of sight and (2) that our line of sight passes through the equatorial plane of the accretion 
flow. 

Our analysis of depolarization is applicable even if the radio emission from Sgr A* is dominated 
by a jet/outflow, rather than the accretion flow as we have assumed. In jet models, it is still natural 
for the highest frequency emission to originate very close to the black hole; in Falcke's model, e.g., 
the J; 100 GHz emission arises from ^ 10 R g , in what is really a "transition region" between the 
accretion flow and the jet (Falcke 1999). In order for this emission to not be depolarized (either in 
situ or propagating through the accretion flow), our constraints on the rotation measure and the 
plasma conditions close to the black hole still apply. 

Two scenarios in which accretion at the Bondi rate could be consistent with observed linear 
polarization at high frequencies are (1) if the high frequency emission arises close to the black hole, 
but in a nearly empty funnel pointed directly towards us (e.g., along the rotation axis of an ADAF) 
or (2) if the high frequency emission from Sgr A* is produced at very large distances from the black 



Bower et al. (1999ab) showed that "bandwidth" depolarization at v ~ 8 GHz requires RM > 10 7 rad m -2 ; this 
constraint does not, however, apply to the depolarization discussed here, namely that due to a spatially varying RM. 
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hole, r J; 10 3 . The former possibility requires a rather special geometry]] and the latter is ruled out 
by the VLBI source size of ~ 10 R g (Krichbaum et al. 1998) and the variability of Sgr A* at ~ 100 
GHz (Tsuboi, Miyazaki, & Tsutsumi 1999). 



4.1. Application to Other Systems 

Although we have have focused our analysis on Sgr A* at the Galactic Center, linear polar- 
ization of high frequency radio emission can be used as a probe of the accretion physics in other 
low-luminosity galactic nuclei (see, e.g., Nagar et al. 2000 for recent VLA observations of LLAGN). 
For a black hole of mass M = mgl{fM & accreting (spherically) at a rate M = 10~ 4 m_4M cc id ~ 
10 23 m_4?ri9 g s , the density, magnetic field strength, and rotation measure in ADAF models are 

n « 3 x 10 6 m_ 4 nig 1 r~ 3 / 2 cm -3 , (15) 

B « 100 mil ™9 1/2 r~ 5 / 4 G, (16) 

and 

RM ns 3 x 10 10 mil m 9 1/2 r ~ ?/4 rad m ~ 2 - ( 17 ) 
Equation ( |i~7| ) shows that large rotation measures and the associated depolarization of synchrotron 
emission by Faraday rotation are generic features of ADAF models (unless m_4 <C 1). 

The absence of observed linear polarization in the radio spectrum of a low-luminosity galactic 
nucleus would be consistent with ADAF models. By contrast, detected linear polarization would 
argue against an ADAF as the source of the observed radio emission. A particularly interesting 
class of systems for future polarimetry are elliptical galaxies in nearby X-ray clusters (e.g., NGC 
4649, 4472, and 4636 in the Virgo cluster). As discussed by, e.g., Fabian & Canizares (1988), Fabian 
& Rees (1995), and Di Matteo et al. (1999, 2000), many of these galaxies have extremely dim nuclei 
given the inferred black hole masses (~ 1O 9 M0) and Bondi accretion rates. Linear polarization 
may shed important light on the physics of these systems. 

We thank Don Backer, Roger Blandford, and Mark Reid for useful correspondence, Bruce 
Draine for useful conversations, and John Bahcall, Heino Falcke, and Feryal Ozel for helpful com- 
ments on the paper. EQ is supported by NASA through Chandra Fellowship PF9-10008, awarded 
by the Chandra X-ray Center, which is operated by the Smithsonian Astrophysical Observatory 
for NASA under contract NAS 8-39073. AG was supported by the W. M. Keck Foundation and 
NSF PHY-9513835. 



7 For example, in Stone et al.'s (1999) simulations of non-radiating accretion flows, the density varies with polar 
angle roughly as p oc sin 2 9 (see also Quataert & Gruzinov 2000); thus for M ~ Meondi, the emission must be confined 
to 6 3° in order for the density to be sufficiently small that high frequency emission is not depolarized. In addition 
our line of sight must lie within ^ 3° of the rotation axis of the flow. 
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A. Faraday rotation in an ultra-relativistic Maxwellian plasma 



Faraday rotation in a cold plasma is described by a change in position angle given by 

dO _ fey ujpUB 



dl 2 u 



(Al) 



where u = ck is the frequency of the radio wave, k\\ is the projection of the wavenumber along 
the magnetic field, u> 2 = Aime 2 /m e is the plasma frequency, and = eB/(m e c) is the cyclotron 
frequency. This corresponds to the usual rotation measure 

6 3 f f rad 

RM = — = — / dl ■ Bn = 2.63 x 10~ 13 x / d\ ■ Bn — =-. (A2) 

A 2 2irm 2 c A J J m 2 y ' 

Here we derive the rotation measure for an ultrarelativistic Maxwellian plasma: 



i?M 7 = 6 \ - A ( dl ■ Bn^J = 2.63 x lO" 13 x / dl ■ Bn^J (A3) 
27rm,gC 4 J 2j z J 2j z m z 

where we have defined 7 = kT e /(m e c 2 ). The dominant correction to the non-relativistic expression 
is the relativistic mass: m e — > 7m e . 

We use the Vlasov equations to calculate the plasma permittivity and hence the dispersion 
relation for electromagnetic waves. For a magnetic field and wavenumber along the z axis, the 
first-order (in the unperturbed magnetic field) permittivity is given by 

m -i4vre 2 eB f j3 1 p\dF m 2 e c 2 



^ 2w m e m e c J (lo — kv z ) 2 p dpp 2 + m 2 c 21 

where the unperturbed distribution function is normalized by J d 3 pF(p) = n, and p\ = p 2 + p 2 . 
For a cold plasma, equation (A4) gives 

Using standard arguments (e.g., Rybicki & Lightman 1979) this leads to the RM for a cold plasma 
given by equation QA2| ) . For an ultra-relativistic plasma, equation ( |A4| ) gives 

where we have not changed the definition of oj p and in the ultra-relativistic regime. Equation 
( |A6| ) for the permittivity gives the ultra-relativistic RM in equation (|A3|). 
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